Singular quasiparticle scattering in the proximity of charge instabilities 



C. Castellani, C. Di Castro, and M. Grilli 

Dipartimento di Fisica, Universita di Roma "La Sapienza", 
Piazzale A. Moro 2, Roma, Italy 00185 



in 

On 



in 



> 
o 



o 
in 



C 

O 

o 



We analyze the behavior of the dynamic scattering ampli- 
tude between Fermi liquid quasiparticles at the Fermi surface 
in the proximity of a charge instability, which may occur in 
the high temperature superconducting cuprates. Within the 
infinitc-U Hubbard-Holstein model in the slave-boson large-N 
technique we find that, in the absence of long-range Coulomb 
forces the scattering amplitude is strongly singular at zero 
momentum transfer close to the phase separation instability 
and it has the same form provided by gauge-field theories. 
In the presence of long-range Coulomb forces the charge in- 
stability occurs at finite wavevectors and concomitantly the 
scattering is still singular but anisotropic. Nevertheless it re- 
mains strong over extended regions of the momentum space. 
In both cases we show how normal state properties are largely 
affected by this scattering. 

PACS:74.72.-h, 71.27.+a, 72.10.-d 



It is generally accepted that the understanding of the 
pairing mechanism in high T c superconductors is related 
to the understanding of the anomalous behaviour of the 
normal phase. 

The anomalous properties of the normal phase have 
been interpreted along two distinct theoretical lines. One 
possible explanation is that the low dimensionality of 
these highly anisotropic systems and their correlated na- 
ture are at the origin of a breakdown of the Fermi liquid 
(FL).In particular the proposal of a Luttinger liquid for- 
mation in two dimension (!]] has been intensively inves- 
tigated ||. The alternative aptitude has been to accept 
the Landau theory of normal FL's as a suitable start- 
ing point. The anomalous properties would then arise as 
a consequence of strong scattering processes at low en- 
ergy between the quasiparticles. Along this line magnetic 
scattering has been considered to be responsible for both 
the anomalous properties of the normal phase and for the 
superconducting pairing || . Strong scattering may even 
lead to a complete disruption of the FL phase. In partic- 
ular it was proposed that excitonic scattering could give 
rise to the so called marginal FL Q , and could also pro- 
vide a pairing mechanism. Singular scattering can also be 
provided by gauge fields |5j , which arise by implementing 
the resonating-valence-bond idea in the t-J model. 

In this letter we want to understand whether phase 
separation (PS) or the incommensurate charge density 
wave (ICDW) instability are sources of strong scatter- 
ing besides the above mechanisms. Indeed the complex 
nature of the phase diagram as a function of doping and 
temperature indicates that various energy scales of differ- 



ent nature (magnetic, excitonic,...) of the same order of 
magnitude compete to determine the low-energy physics 
and may lead to various instabilities, among which PS or 
charge instabilities may play a relevant role. 

After PS was shown to be present in the phase di- 
agram of the t-J model we pointed out that PS 
commonly occurs in models with short range interaction 
|^-|l5|, provided the strong local e-e repulsion inhibits 
the stabilizing role of the kinetic energy. We therefore 
stressed that PS and superconductivity can be related 
phenomena irrespective of the nature of the short-range 
interaction. Emery and Kivelson [|l6| suggested that, al- 
though long range Coulomb (LRC) forces spoil PS as a 
static thermodynamic phenomenon, the frustrated ten- 
dency towards phase separation may still be important 
and give rise to large amplitude collective density fluctu- 
ations. Approaching the problem within a coarse-grained 
model, they suggested that these fluctuations may be 
responsible for the anomalous behaviour of the normal 
phase and for the superconducting pairing. 

To assess the relevance of PS as a mechanism for 
anomalous scattering, we here determine the dynami- 
cal effective scattering interactions among quasiparticles 
close to a charge instability, both in the presence and in 
the absence of LRC forces. We carry out this analysis 
within a microscopic treatment of the Hubbard-Holstein 
model in the infinite- U limit. We find that, both in 
the presence and in the absence of LRC forces, the dy- 
namic effective interaction turns out to have a singular 
behaviour, strongly affecting the single-particle and the 
transport scattering time. 

-The model - Although our results are quite generic 
of models with PS, to be specific we use as a simple 
paradigm the two-dimensional Hubbard model with an 
additional dispersionless phonon mode A coupled a la 
Holstein 
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where and indicate nearest-neighbor and 

next-nearest-neighbor sites respectively and J2a = 
J2 a c ia c iv i s the local electron density, which in mo- 
mentum space is given by p q = J2k,a c l+ q ,a c k,v- Tne 
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last term of Eq.(l) is a Coulombic potential between 
electrons on a two-dimensional square lattice (with lat- 
tice spacing a in the x and y directions) embedded 
in a three-dimensional space with a separation d be- 
tween the planes in the z direction. The dielectric con- 
stants in the plane and perpendicular to it are en and 
e_L respectively and the Coulombic coupling constant 
Vc = e 2 d/2e±a 2 . The momentum dependence of the 
potential on the z=0 plane is found to be G(q x ,q y ) = 
e±(a/d) 2 ( cos ( a( fa) + cos(aq y ) — 2) — 1. As usual, the sum 
in the Coulombic potential does not include the zero- 
momentum component, since we are supposing that the 
diverging q = interaction between the electrons is can- 
celed by the contribution of a uniform positively charged 
ionic background. 

Since we are interested in strong local repulsion we take 
the limit U — > oo, which gives rise to the local constraint 
of no double occupancy n ^ < L To implement this 
constraint we use a standard slave-boson technique , 
by performing the usual substitution c| CT — > c\ a bi, Ci a — > 
b\ci a and introducing a Lagrange multiplier field A*. The 
quartic Coulombic term can be decoupled by a Hubbard- 
Stratonovich transformation introducing an additional 
real bosonic field Yi. Within the large- TV expansion we 
assume that the spin index runs from 1 to N and we re- 
lax the constraint to the form ^2 a c ia .Ci a + &|i>j = •y. A 
suitable rescaling of the hoppings t — > t/N and t' — > t' /N 
must, in this model, be joined by the similar rescaling of 
the e-ph coupling g — > g/^/N and of the Coulomb in- 
teraction Vc — > Vc/N in order to compensate for the 
presence of N fermionic degrees of freedom. 

The model can first be solved in the mean field (N = 
oo) approximation by setting the 6,j and Aj bosons to 
their constant self-consistent values b and A respec- 
tively. The system then results into free quasiparticles 
with a shifted chemical potential /i = fiQ — Ao and a dis- 
persion Ef. — — 2ir§£fc with Ek = (cos(afe x ) + cos(ak y )) + 
(f /t) (cos(ak x + ak y ) + cos(ak x — ak y )) where r 2 . = 
b^/N = 8/2. For any finite doping 5 at T=0 the sys- 
tem is a Fermi-liquid where the mean-field value of the 
slave-boson field bo multiplicatively reduces the hoppings, 
t — ► f&Q, t' — > t'b 2 ,, thus enhancing the effective mass of 
the quasiparticles. At this level the mean-field self-energy 
does not introduce a finite quasiparticle lifetime. 

The effective interaction leading to scattering between 
quasiparticles arises from the exchange of the bosonic 
fields in the 1/N corrections beyond mean field approx- 
imation. One can define a four-component field — 
(Sr, SX, <p, Y). <j>i = {A\+Ai)/{2VN) is the lattice dis- 
placement field, and 5ri and SXi are the fluctuating part 
of the &i-field amplitude and of the the Lagrange mul- 
tiplier respectively. The leading-order expressions of the 
effective scattering amplitude in the particle-hole channel 
can be written as 



T(k, k'; q,uj) = -Y,^ (*', -q) D» v (q, lo) A" (k, q) . (2) 

where A'' are the vertices coupling the fermionic quasi- 
particles to the bosons, A r (k,q) = 
-2tr 2 (e k+q/2 +Ek- q /2), A A (fc,g) = i, ^<f>(k,q) = -2g 
and A Y {k,q) = i. D> iv (q,uj) = (a^(q, u)a v {-q, -u)) = 
N^ 1 (2B + H(q,uj))~l is the leading order boson propa- 
gator with self-energy corrections given by the fermionic 
bubbles (which include quasiparticle-boson vertices) . 
B is the bare boson-propagator matrix. An expression 
similar to Eq. (||) holds for the scattering amplitude in 
the particle-particle channel. 

- The results - The evaluation of the density correlation 
function P(q,u) = {1/N) w)n ff /(-g, -u)) 

provides information on the stability of the system. In 
particular a divergence in the static density-density cor- 
relation function P(q, lo = 0) signals the occurrence of 
PS (at q — > 0) or of CDW instabilities (at finite q's). A 
complete investigation of the static and dynamical prop- 
erties of the present model together with the analysis of 
its stability was already carried out in a previous work 
|l5| . Here we just mention that, within the present for- 
malism, the model ([!]) displays a phonon-driven charge 
instability even for rather small e-ph coupling. In the 
absence of LRC forces this instability occurs before any 
other finite q instability up to intermediate-large doping. 
At large doping the instability requires larger e-ph cou- 
plings and occurs at finite q w 2kp in some directions 
[in particular (1,0) and (0,1)] signalling the occurrence of 
incommensurate CDW. The introduction of LRC forces 
eliminates the small-g divergence in the static correlation 
function always giving rise to finite q instabilities. The 
critical q, in this case, is not related to any pseudonest- 
ing of the Fermi surface but it depends on the strength 
of the LRC forces and on the momentum dependence of 
the poles in the divergent static correlation function with 
only short-range forces. Pair formation is always found 
near the instabilities. 

As already announced, we want to report on the be- 
havior of the quasiparticle scattering close to both the 
PS and the CDW instability. It is worth noting that a 
divergent scattering amplitude will follow from a diver- 
gent correlation function P. Indeed a divergent boson 
propagator enters in the expression of both quantities 
establishing a clear connection between the charge insta- 
bility and the singular quasiparticle scattering. 

In this regard we carried out an extensive analysis of 
the real and imaginary parts of the scattering amplitude 
between the quasiparticles on the Fermi surface. 

As expected, near the PS instability (Vc = 0), the 
anomalous behavior of T is identified to be of the form 

T(q,ui) « T < (3) 

K ' Bq 2 -iuj£+D W 

As shown in Fig. la, for the model with t = 0.5eV, 
t'/t = -1/6, ujq = 0.04eV and g = 0.194eV, D = D(S — 
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S c ) vanishes linearly when, for a given g, the instability 
takes place at the critical doping S c = S c (g). For the 
same values of the parameters, Fig. lb reports the real 
static scattering amplitude r(q,w = 0) as a function of 
transferred momentum q at various dopings close to the 
q = instability. In this case the singular part of the 
scattering amplitude displays a quite isotropic behavior 
(at q c = 0). 

In the absence of LRC forces the singular behavior of 
T(q — > 0, u) = 0) at the PS instability is by no means sur- 
prising within a FL framework. Indeed, the FL expres- 
sion for the compressibility is n = 2v* / (1 + 2v*r uJ ) = 
2f*(l — 2is*T q ), where T w and T q are the standard dy- 
namic (u> — > 0, q — 0) and static (q — > 0, u> = 0) limit of 
the scattering amplitude. This indicates that a divergent 
k, when the quasiparticle mass remains finite (y* < oo), 
only happens when Fq = 2v*T UJ — * — 1 (Pomeranchuk 
criterion). At the same time if r g — * oo. We like to point 
out here that the above arguments keep their full validity 
irrespective of the mechanism leading to PS. 

In the presence of LRC forces the singular part of T 
can be written as 



(4) 



where toq = D' + B'\q — q c \ 2 The behavior of the 
mass term D 1 is shown in Fig. 2a for the model with 
t = 0.5eV, t'/t = -1/6, V c = 1-lcV, uj q = 0.04eV 
and g = 0.240eV as a function of S — 5 C . For these pa- 
rameters the instability first occurs at 5 C = 0.194 and 
q c w (±0.28/a,±0.86/o), or q c w (±0.86/a, ±0.28/a). 
Analogously to Fig. lb, Fig. 2b displays the strong dop- 
ing dependence of the static scattering amplitude as a 
function of momenta in the q c direction. However, as 
shown in Fig. 3, the scattering is quite strong, although 
non-singular, in all directions for \q\ w |q c |. We also 
checked that the (almost) isotropic contribution to the 
static scattering amplitude it is much less fragile under 
doping variations. 

The imaginary term in the denominators in the r.h.s. 
of Eqs.(||) and (jj) reproduces on a wide range of 
transferred momenta q the behaviour of the imaginary 
part of the mean field fermionic polarization bubble 
Im \P°(q, w)] oc ui/q at small uj. This indicates, that, 
despite the complicated formal structure of the scattering 
amplitude (|J) arising from the matrix form of the boson 
propagator , near the instability a simple RPA-like 
structure results in the final expression. The forms (|^) 
and (^ arc, therefore, generic of PS or charge instabili- 
ties. 

It is apparent that the effective interaction (||) in the 
isotropic pure-short-range case, has the same form as the 
one mediated by gauge fields if it were not for the 
doping-dependent mass term and for the completely dif- 
ferent origin of the singularity. The LRC forces make the 
scattering anisotropic rendering it more similar to the 



magnetic fluctuation case , even though the really sin- 
gular behaviour is for differently oriented momenta and 
large scattering is present in any direction. 

Near PS at 5 — S c , the electrons are coupled to a mode, 
which is soft at long wavelength, and the same results of 
Rcf. Q apply. Therefore in two dimensions the inverse 
scattering time is proportional to T 2 / 3 and the resistiv- 
ity p is proportional to T 4 / 3 with a crossover temperature 
T* to the standard FL behavior, which in our case turns 
out to be T* oc (6 — S c ) 3 ^ 2 . The resistivity is anomalous 
but it is not linear in T as in the optimally doped high 
T c materials. An additional the limitation of the pure 
short-range case may be the difficulty in approaching the 
instability line which might be embedded in a region sta- 
bilized by the Maxwell construction. 

In the more physical case including LRC forces, the 
proposal of associating the anomalous properties of the 
normal phase to the presence of a T = quantum critical 
point [ fj"9|]20| has a natural realization within the charge- 
instability mechanism. The outcome for the transport 
properties should be similar to the one described within 
the nearly antiferromagnetic FL scenario H, i.e. a linear 
resistivity for temperatures larger than a crossover tem- 
perature T* , which, in our case would be proportional 
to D' in Eq.(j4|). However, this result was questioned in 
Ref. [^l| for the AF fluctuations because only few "hot" 
points on the Fermi surface are connected by singular 
interactions, with an anomalous inverse scattering time 
proportional to VT. Generically the resistivity would be 
dominated by an inverse scattering time proportional to 
T 2 , then linear resistivity would only appear for T>r. 
In our case, instead, a strong rather isotropic scattering 
persists rendering the above limitation less effective. 

The above scenario connecting charge instabilities to 
the anomalous normal properties of the superconducting 
cuprates raises the problem of the relevance of charge 
fluctuations in these systems j22|. Indeed commensurate 
CDW have been recently observed (23) in a related com- 
pound (La2- a; _2 / Nd2 / Sr 2 ;Cu04 with x = 0.12 and y = 0.4) 
possibly due to the pinning of electronic dynamical CDW 
by a suitable underlying lattice structure. It is therefore 
quite important to assess the electronic origin of this 
CDW, e.g. by detecting non-linear effects like electric- 
field-induced depinning. This would strongly support the 
existence of dynamical CDW fluctuations when pinning 
is not effective leaving the materials in the metallic phase. 
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FIG. 1. (a) Mass D as a function 8 - S c for t = 0.5eV 
t' = -1/6*, Vc = 0, loo = 0.04cV and g = 0.194eV. (b) Static 
scattering amplitude for the same parameters as in (a) as a 
function of the transferred momentum q in the (1,0) direction. 
The doping 8 = 0.2, 0.205, 0.22, 0.25, 0.28 increases from the 
lower solid line to the upper solid line. 
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FIG. 2. (a) Mass D' as a function 8 - 8 C for t = 0.5eV 
t' = -1/6*, Vc = 1-leV, ujo = 0.04eV and g = 0.312eV. 
(b) Static scattering amplitude for the same parameters 
as in (a) as a function of the transferred momentum q 
in the q c w (±0.28/a, ±0.86/a), direction. The doping 
8 — 0.195,0.2,0.205,0.22 increases from the lower solid line 
to the upper solid line. 
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FIG. 3. Momentum dependence of the static scattering am- 
plitude for the same parameters as in Fig. 2a at S = 0.195. 
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